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Abstract 

In this paper, we extend the analyses carried out in a previous article for VKVF-fusion 
to the case of Higgs production via ZZ-fusion within the Standard Model at the Next 
Linear Collider, in presence of electromagnetic radiation due real photon emission. Cal- 
culations are carried out at tree-level and rates of the leading order (LO) processes 
e + e~ — ► e + e~H — ► e + e~bb and e + e~ — ► e + e~H — ► e + e~WW — ► e + e~jjjj are compared 
to those of the next-to-leading order (NLO) reactions e + e~ — * e + e~H{^) — > e + e~bb"f and 
e + e~ — > e + e~H( r y) — > e + e~WW( , y) — ► e + e~jjjj7, in the case of energetic and isolated 
photons. 
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In a previous paper |j] we studied the effects of hard photon emission in the V^W-fusion 
reactions 

e + e" -> v e v e H{i) -> u e v e blrf, (1) 
e + e~ -> v e v e E{^) -> p e u e WW(j) -> P e ^ejjjj7- (2) 

The production mechanism e + e~ — ► v e v e H and the decay channels U — > 66, jjjj represent some 
of the most likely signatures of Higgs processes at the Next Linear Collider (NLC) @-||, within 
the Standard Model (SAi), for a H scalar in the intermediate (1) and heavy (2) mass range 
(IMR and HMR, respectively )Q. 

In that study we focused our attention to the case of hard and isolated photons. The 
motivation of that analysis came from the observation that the role of energetic photons from 
the initial state is expected to be determinant at the e + e~ colliders of the next generation. 
In particular, whereas at LEP1, SLC and partially LEP2, the finite width of the gauge boson 
resonances heavily suppresses events with hard photons produced in the e + e~ annihilation 
subprocess, at a high energy NLC (when y/s ~ 300 GeV) this is no longer true. Furthermore, 
since the energy of the beams is much larger, the probability that the incoming electrons and 
positrons radiate is somewhat bigger. Therefore, one naturally expects the NLC to produce 
rather copiously events accompanied by hard electromagnetic (EM) emission and these cannot 
be ignored in realistic phenomenological analyses. 

In many instances, one can give account of the Initial State Radiation (i.e., bremsstrahlung of 
photons from the incoming electron/positron lines) by means of the so-called 'electron structure 
functions' ||. In addition, for certain designs of the NLC, the effects due to such a radiation are 
dominant with respect to those due to beamsstrahlung and Linac energy spread phenomena 
0: such that, in phenomenological studies, one can consistently deal with bremsstrahlung 
photons only. However, on the one hand, this approach is realistically applicable only in the 
case of annihilation- and conversion-type Feynman diagrams: that is, when the e ± -lines are 



connected to each other via s- or t, u-channel interactions p0|| . When the electron and positron 
lines are disconnected by gauge boson currents, like in the case of lyVF-fusion processes, a 
separation of the radiation emitted by the incoming fermion lines from that radiated by the 
virtual VF-boson lines (see Fig. 1-2 in Ref. 0) is indeed not gauge invariant. On the other 
hand, in experimental data, radiative events with photons coming from the Higgs production 
mechanism are not distinguishable from those in which they come from the Higgs decays. 
To give reliable predictions for the corrections at the order 0(a em ) to the leading processes 



e e 



v e v e E -> v e v e bb, (3) 



2 For a NLC with y/s ~ 500 GeV and Mh ~ 2M\y the production rates of the bremsstrahlung process 
e + e~ — > ZH are comparable or larger than those of the VKW^-fusion channel || u\. 
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e+e -> Z/ e Z/ e ^ -> f^e^W -> ^ejjjj, (4) 

one should compute the contributions due to real photons as well as those due to virtual photons 
in the loop diagrams. Their summation through the orders 0(al m ) and 0(a 7 cm ) would allow one 
to cancel the infrared (soft and collinear) divergences appearing in the matrix elements (MEs). 
However, to compare the tree-level rates of the LO processes (3)-(4) with those of the NLO 
ones (l)-(2) for hard photons only, allows one to assess whether (at higher order) complications 
should be expected in working out the position of the resonances and/or in establishing their 
line-shape, or whether the overall effect of the complete 0(a em ) corrections will be in the end 
only matter of a different normalisation (of the LO spectra). 

In the approach of Ref. [jlj, virtual photon contributions were ignored, for both the cases 
of what we called 'production radiation' and 'decay radiation'. (In detail, the first is due to 
photons emitted during the Higgs production mechanism, whereas the second refers to photons 
radiated in the Higgs decays, and they are separately gauge- invariant.) As loop diagrams were 
missing, the mentioned summation was not performed there. In that paper, we compared the 
rates of the lowest order processes (|3|) and (f|) to those at next-to-lowest order from reactions ([]]) 
and (0) over the phase space regions defined by the cuts: j ' 7 > 1 GeV and cos 6^ - n < 0.95. 
The aim was to make the point that in many instances (that is for various combinations of 
NLC centre-of-mass (CM) energies, a/s, and Higgs masses, Mh) the rates of the hard radiative 
processes can be rather large compared to those of the leading reactions and that in some cases 
a complete C(a em ) calculation would be desirable in order to assess the relative importance of 
events of the type (§) with hard photons, with respect to those in which the electromagnetic 
emission is either infrared or virtual and to the lowest order rates as well, as a shift in the position 
of the Higgs peaks could be expected at NLO (see again Ref. Q for a fuller discussion). Another 
feature of interest to phenomenological analyses stressed in Ref. [0J is the fact that rather clear 
Higgs signals can appear in the spectra of the invariant masses M 6 g 7 and Mjjjj-y- The last point 
outlined in Ref. JE| is that it is generally impossible to separate in radiative events dH)-(D the 
'production radiation' from the 'decay radiation', at least by exploiting the spectra in photon 
energy, E^, and in missing transverse momentum of the neutrino pair, p™ lss . 

We intend in this brief report to perform studies similar to those of Ref. ]I] for the case 
of a SM. Higgs boson produced via the ZZ-fusion mechanism and decaying through the two 



mentioned channels, by studying the reactions: 

e + e _ — > e + e~H(j) — > e + e~6&7, (5) 

e + e~ — > e + e~H(~f) — > e + e~WW(-y) — > e + e"jjjj7, (6) 

e + e~ — > e + e~ H — ► e + e~bb, (7) 

e + e~ -»• e + e~H -> e + e~WW -> e + e~jjjj. (8) 
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As in this case infrared divergences can occur also in the bremsstrahlung of soft and collinear 
photons off the e ± -lines in the final state, we impose that cos 9 e /tr/,e/ j-y < 0.95 (that is, separation 
of the photon also from the outgoing electron/positron). Moreover, since the e ± -particles in 
the final states are detectable (contrary to the neutrinos in the VFW-fusion) , we ask p^ b ' ^ > 1 
GeV. The Feynman diagrams describing the reactions (^)-(|§) at tree-level are the same pictured 
in Fig. 1-2 of Ref. f]J, with graphs 5-6 there replaced by those in Figs, la-b of the present 
paper. Their calculation has been carried out in the same way as in Ref. |l|], so we refer to that 
publication for all details. 

We have verified that also in the case of the ZZ-fusion processes the spectra in photon 
energy and in missing transverse momentum of the e + e~ pair show a similar behaviour to that 
described in case of processes (|I])-(fD. Therefore, we concentrate here only on the case of the 
integrated rates and in the invariant mass spectra M b i, M^, M fe g 7 and Mjjjj 7 . Our results are 
presented in Figs. 2a-b and Figs. 3a-b. For consistency, the choice of the collider CM energies 
and Higgs masses we have adopted here is the same as in Ref. |l|. 

Many of the features typical of processes ([TJ)-(|]) are reproduced in case of processes ([S])- 
(|8|) as well: yet, differential and integrated rates are here smaller by an order of magnitude 
at least (compare Figs. 2a-b to Figs. 3-4 of Ref. IJ), both at LO (as well known) and at 
NLO. Concerning the total cross sections, one can notice (Fig. 2a) that in the IMR the NLO 
contribution due to hard photons is never greater than « 10% of the LO rates, for all yfs = 
300, 500 and 1000 GeV, when Mh varies in the IMR. The same can be affirmed for a scalar 
with mass in the whole of the HMR (Fig. 2b). 

Figs. 3a-b show the distribution in the invariant masses of the Higgs decay products: i.e., 
M fe 5( 7 ) and M^^y Things are here very similar to the case of the corresponding WW- fusion 
processes, such that the same comments contained in Ref. [l| hold here too. However, a 
different feature of ZZ-fusion processes compared to the HT^-fusion ones is that the relative 
contribution of graphs 5 and 6 to the total cross section at NLO is smaller in the former. This 
is due to the additional cut between the final state electron/positron and photon directions 
implemented in case of processes (fjH(|J), as the EM emission along the e ± lines tends to be 
highly collinear. This ultimately reduces the relevance of NLO hard photon events, compared 
to the lowest order rates. In fact, for the bb invariant masses (left plots in Fig. 3a), the smearing 
towards low values (that is, M h i < M#) of the Higgs resonances due to the hard EM radiation 
in the H — > bbj three-body decays induces an effect of only « 1% at the most, for all Higgs 
masses and collider energies. The suppression of the 'decay radiation' is once again visible also 
in the invariant mass of the 667 system (plots on the right in Fig. 3a), because of the absence 
of clear Breit-Wigner peaks |TJ. For a Higgs scalar in the HMR, we notice that the shape of 
the mass distributions on the left hand side of Fig. 3b is rather different at the two orders. The 
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effect is visible at all energies: for example, see the open and shaded dotted histograms (see also 
Ref. 0). However, the quantitative effect is rather small compared to the case of lyVF-fusion. 
Finally, Higgs peaks clearly appear in the Mjjjj 7 spectra (right hand side of Fig. 3b). 

Therefore, we conclude by remarking that the corrections due to hard photons in reactions 
induced by ZZ-fusion are generally smaller compared to the case of VFW-fusion processes, 
both in the intermediate and heavy mass range of M#. In particular then, when folding the 
real radiation computed here with the virtual one and with the lowest order rates to produce 
the complete result at the 0(a em ) higher order, the smearing of the resonant Higgs peaks 
should be less relevant than in the W^VF-fusion case, and negligible in general. However, in the 
HMR and especially for very large NLC energies (i.e., around and above 1000 GeV, see Fig. 3b 
bottom left), a complete NLO calculation would be desirable for phenomenological analyses 
that consider ZZ-fusion on its own. We finally notice that for ZZ-fusion Higgs processes (like 
for H^H^-fusion) the inclusion of four-jet decays via the ZZ-channel (i.e., H — > ZZ — ► j j j j (t) ) 
should reduce the relevance of hard radiative contributions; and that a careful treatment of 
the Higgs-strahlung process e + e~ — > ZH (followed by Z — > e + e~) and of its interference with 
the ZZ-fusion mechanism would be needed (see the discussions in the last two paragraphs of 
Section 3 in Ref. |l[). The last aspect is especially important when dealing with processes 
(5)-(8), as the LO rates of the two Higgs mechanisms (i.e., Zff-production and ZZ-fusion) are 
comparable for all combinations of y/s and Mg considered here (whereas VFVF-fusion always 
dominates for yfs ~ 500 GeV). 
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Figure Captions 

Fig. 1 New Feynman diagrams at tree-level which substitute those corresponding to the WW- 
fusion channel when the photon is emitted by W bosons in the Higgs production mecha- 
nism (see diagrams 5-6 in Figs. 1-2 of Ref. ]T[): (a) in case of process (|5J); (b) in case of 
process (||). 

Fig. 2 Cross sections of the processes: (a) @ and (0); (b) (|) and (§); as a function of the 
Higgs mass, at y/s = 300 GeV (upper plot), -^/s = 500 GeV (central plot), yfs = 1000 
GeV (lower plot), after the cuts pj! 1 ' 1 > 1 GeV and cos^//^ < 0.95, where h = b or j, as 
appropriate. 
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3 Distributions in: (a) invariant mass of the bb (left plots) and of the 667 (right plots) 
systems for processes <g) and (0), for M H = 60,100,140 GeV at ^ = 300(500) [1000] 
GeV; (b) invariant mass of the jjjj (left plots) and of the jjjj'y (right plots) systems 
for processes (§) and ©, for M H = 180, 220, 260(180, 260, 340) [180, 300, 420] GeV at 
= 300(500) [1000] GeV. The following cuts have been applied: > 1 GeV and 

cos# e /h 7 < 0.95, where h = b or j, as appropriate. The values of the CM energy are: 
y/s = 300 GeV (upper plots, bins of 2 GeV); y/s = 500 GeV (central plots, bins of 4 
GeV); y/s = 1000 GeV (lower plots, bins of 5 GeV). In the left plots the upper histograms 
refer to rates from the non-radiative processes (^) and (§), whereas the lower histograms 
correspond to rates from the radiative processes (|5]) and @. The latter are shaded. 
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Diagrams by MadGraph 




Fig. 1a 



New graph 5 New graph 6 




Fig. 1b 



New graph 5 New graph 6 
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e + e -» e + e jjjj {7) 
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p T 7 > 1 GeV 



e e -> e e bb (7) 
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, p T > 1 GeV 



e + e -> e + e jjjj (?) 
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